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Abstract 

We present a phase diagram of the different kinds of congested traffic that 
are triggered by disturbances when passing ramps or other spatial inhomo- 
geneities of a freeway. The simulation results obtained by the nonlocal, gas- 
kinetic-based traffic model are in good agreement with empirical findings. 
They allow to understand the observed transitions between free and various 
kinds of congested traffic, among them localized clusters, stop-and-go waves, 
and different types of "synchronized" traffic. We also give analytical condi- 
tions for the existence of these states which suggest that the phase diagram is 

universal for a class of different microscopic and macroscopic traffic models. 
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For some years now, the various observed nonlinear states of freeway traffic flow, in- 
cluding traffic jams and stop-and-go waves, have attracted the interest of a rapidly growing 
community of physicists [0,0- The same applies to simulation approaches that allow to re- 
produce the transitions between the different collective states of motion . A particular 
attention received the recent discovery of a first-order phase transition to "synchronized" 
congested traffic |||7| , which stimulated an intense investigation of the complex phenomena 
associated with ramps, yielding an explanation of the hysteretic phase transition |||5|]. It is 
also known from other effectively one-dimensional many-body systems with short-range in- 
teractions that localized inhomogeneities can cause a variety of phenomena, including phase 



transitions M and spontaneous symmetry breaking [10 



To characterize the parameter-dependence of the possible states of a system resulting in 
the long run, phase diagrams are a very powerful method. They are of great importance 
in thermodynamics with various applications in metallurgy, chemistry, etc. Moreover, they 
allow to compare very different kinds of systems like equilibrium and nonequilibrium ones, or 
microscopic and macroscopic ones, whose equivalence cannot simply be shown by transfor- 
mation to normal forms |TT|JT2]| . Defining universality classes by mathematically equivalent 
phase diagrams, one can even classify so different systems as physical, chemical, biological 
and social ones, as done in systems theory. 

In the following, we will present and explain the phase diagram that we obtained for 
the nonlocal, gas- kinetic-based traffic model ||, when we studied a freeway with a ramp 
in the presence of a single perturbation [Fig. [p. By systematical variation of the inflow 
at the upstream freeway boundary and the ramp, we were able to reproduce the observed 
transitions between different traffic states, some of which have not been explained before. 
We will also present analytical conditions for the existence of the various traffic states. Since 
we are interested in the generic properties of the model, we will formulate it in dimensionless 
quantities jn§ by measuring times in units of the relaxation time r « 40 s, and distances 
in units of the inverse l/p max of the maximum vehicle density p max ~ 140 vehicles/km. The 
equation for the vehicle density p(x, t) at position x and time t is a continuity equation with 
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a sink/source term and reads 



dp d{pV) = grmpQM) 

dt dx L rmp 



V denotes the average vehicle velocity, L Tmp is the length of the ramp, Q rmp is the flow 
of vehicles entering the freeway (Q rmp > 0) or leaving it (Q rmp < 0) along the ramp of 
length L r mp> divided by the number n of freeway lanes. The velocity equation contains a 
convection term (due to the movement of the velocity profile with velocity V), a pressure 
term (reflecting dispersion effects due to a finite velocity variance 6), a relaxation term 
(describing an adaptation to a desired velocity Vo), and an interaction term (corresponding 
to braking maneuvers) [|13|]: 

ev +v av = _i_ a j^ 

dt ox p ox 



-PA(P,)^%B(5 V ). (2) 



Vo has a meaning similar to the Reynolds number (since it allows for traffic instabilities 
when it exceeds the value V c = 61). In addition, P is a scaled cross section. 

A(p) = 0.171 + 0.417{tanh[10(p - 0.27)] + 1} (3) 

represents a structure factor, which is of order unity and determines the variance via the 
constitutive relation 9 = A(p)V 2 , but also the form of the flow-density relation in equilibrium 
@ (see Fig. |). Finally, 



B(S V ) 



5 v ^= + (l + 5 2 v ) dy 



2tt J-oo \/1tx 



(4) 



with the dimensionless velocity difference by = (V — V a ) j\[B + 6 a is a Boltzmann factor 
arising from the vehicle interactions. An index "a" indicates that the respective quantity is 
evaluated at the advanced "interaction point" x a = 37 + 7(1 +TV) rather than at the actual 
position x. The related nonlocality reflects the anticipative driver behavior and is essential 
for the realistic properties of the model and its robust and efficient numerical solution. 7 is 
an anticipation factor and T about the safe time headway. In our simulations, we used the 
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parameter values Vo = 171, P = 0.31, 7 = 1.2, and T = 0.043, which are typical for Dutch 
freeways ||13|| , but our results are not very sensitive to their particular choice. 

The nonlocal, gas-kinetic-based traffic model reproduces the characteristic properties of 
traffic flows formulated by Kerner and Rehborn [rj|. Furthermore, for a homogeneous 
freeway without ramps, it shows the phase ( "instability" ) diagram that was postulated 
by Kerner and Konhauser on the basis of macroscopic simulations For a given density p, 
there exists a state of homogeneous traffic with equilibrium velocity V e (p) and equilibrium 
flow Q e (p) = pV e (p) (see Fig. H), and there are four critical densities p ci . For densities 
p < pd and p > p C 4, homogeneous traffic is stable with respect to localized perturbations, 
and for a range p C 2 < p < p c % of intermediate densities, it is linearly unstable, giving rise to 
cascades of traffic jams ( "stop-and-go traffic"). For the two density regimes p c i < p < p C 2 
and p C 3 < p < p C 4 between the stable and the linearly unstable regions, it is metastable, i.e., 
it behaves nonlinearly unstable with respect to perturbations exceeding a certain critical 
amplitude, but otherwise stable. For the discussion of the phase diagram, an additional 
range p cv < p < p c % of convective stability will be important, where homogeneous traffic is 
linearly unstable, but the growing perturbations are convected away from any fixed location 
0. 

We will denote by Q ci = Q e (p C i) the equilibrium flows corresponding to the critical 
densities introduced above. Moreover, there exists a characteristic outflow Q Q ut from traffic 



jams, stop-and-go waves, etc. [13 , 14| , that is nearly independent of the surrounding density 
and the type of congested traffic, at least if there is no ramp. If traffic relaxes at the location 
of an on-ramp, the observed outflow is Q ut(Qrm P , i mp ) < Qout, but we find Q out — > Qout 
for L rmp — ► 00 or Q Tmp — > 0. For our model parameters, we obtain p cl = 0.11, p c2 = 0.12, 
p cv = 0.34, p c3 = 0.36, and p C 4 = 0.40, related to Q c i = 16.2, Q c2 = 17.5, Q cv = 14.6, 
Q c3 = 14.1, and Q c4 = 13.2. If L > 50, we obtain Qout ~ Qc2, otherwise Q out decreases with 
growing Q rmp and decreasing L rmp . For these values (which, with the exception of Q ou t, 
are determined from situations without ramps!), the equations for the phase boundaries 
suggested below agree quantitatively with the ones in the phase diagram in Fig. [I]. 
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In our simulations, we investigated the congested traffic states that formed near an on- 
ramp when passed by a perturbation on the freeway (Fig. 0). For the perturbation we chose 
a fully developed density cluster that did not change its amplitude or shape anymore and 
travelled upstream with constant speed. When the density cluster reached the ramp, it 
induced different kinds of congested states (Figs. [I] through H), depending on the inflows Q in 
and Qrmp at the upstream boundary of the freeway and the ramp, respectively. 

Now, we will explain the different states in the phase diagram starting with relatively 
high ramp flows. When traffic breaks down (a sufficient criterium for this is given by 
(Qout + Qrmp) > Qmax, where Q max denotes the maximum of the equilibrium flow), we 
find the formation of a growing, extended region of congested traffic with average flow 
Qcong = (Qout — Qrmp) in front of the downstream end of the on-ramp ||. Because of the 
fixed downstream end, any perturbation will be convected out of the congested region, if 
the homogeneous solution of density p con g and flow Q con g = Qc(Pcong) is convectively stable, 
i.e. if Qcong < Qcv In this case, we have homogeneous congested traffic (HCT) [Fig. 0(a)], 
which corresponds to an equilibrium solution on the high-density branch of the flow-density 
diagram [Fig. 0(a)]. Otherwise, oscillatory congested traffic (OCT) forms [Figs. |||(b) and 
H(b)]. Thus, the boundary between HCT and OCT is given by 

HCT-OCT: Q rmp = Qout - Qcv (5) 

The oscillation amplitudes grow with decreasing ramp flow, until they reach the low-density 
part of the flow-density diagram associated with free traffic [Fig. @(c)]. This means that we 
have an alternation between free and congested traffic which defines the state of triggered 
stop-and-go waves (TSG) [Fig. 0(c)]. It turns out, that TSG is characterized by Qtsg > 
Qtsgj where Qtsg denotes the average flow of the TSG state. Hence, the boundary between 
TSG and OCT is given by 

OCT-TSG: Q rm p = Qout — Qtsg • (6) 
However, the precise value of Qtsg * s hard to determine. 
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The minimal inflows to sustain the TSG states can be derived from the triggering mech- 
anism illustrated in Fig. ^(c). When passing the ramp, the initial localized perturbation 
causes a secondary perturbation travelling downstream |T5|. Since the amplitude of this 



secondary perturbation is always small, regardless of the amplitude of the primary pertur- 
bation, it only grows on the instability condition Qdown > Qc2, where Qdown = (Qm + Qrmp) 
denotes the flow downstream of the ramp. With growing amplitude, the triggered perturba- 
tion changes its propagation speed, reverses its direction, and finally induces another small 
perturbation when passing the ramp, etc. In this way, the interplay of the perturbations with 
the ramp defines an intrinsic timescale and wavelength of the triggered stop-and-go waves. 
If the above instability condition is not fulfilled, instead of TSG we find a single moving 
localized cluster (MLC) [Figs. |3|(d) and 0(d)], similar to the transition between stop-and-go 
waves and localized clusters in homogeneous traffic ||]. The respective boundary is given by 

TSG-MLC: Q rmp = Q c2 - Q m . (7) 

The MLC continues the course of the initial perturbation which, however, changes its propa- 
gation speed at the ramp, probably because of Q out < Q on t- If Qm < Qci, the upstream flow 
cannot sustain localized clusters, since it is stable to perturbations, and the MLC becomes 
a pinned (standing) localized cluster (PLC [or SLC]) [l5j at the on-ramp [Fig. 0(e)]: 

MLC-PLC: Q in = Q cl . (8) 

A pinned localized cluster will grow and, thereby, give rise to a an extended form of 
congested traffic (CT), if the sum (Q in + Q Tmp ) of the inflows exceeds the self-organized 
outflow Qout from CT behind ramps. The corresponding boundaries are given by 

PLC-CT: Qrmp = Qout — Qin • (9) 

They are the same for HCT and OCT, so that there could be also stationary and oscillating 
variants of PLC (SPLC and OPLC). 
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Finally, no congested state can be formed, if Qdown < Qci, since traffic flow is stable 
to any perturbation, then. This defines the boundary between the pinned localized cluster 
state and free traffic (FT): 

PLC-FT: Q rmp = Q cl -Q in . (10) 

Due to the metastability between p c \ and p c2 , the transition to free traffic depends on 
the amplitude of the perturbation. For small perturbation amplitudes, stable traffic flow 
associated with free traffic expands up to Qdown < Qc2- Hence, the MLC and PLC regimes 
may disappear in cases of small disturbances, leading to the boundaries FT-TSG, FT-OCT, 
and FT-HCT. In conclusion, the transition between free and congested traffic is of first order 
(i.e. hysteretic) 0, while the other transitions seem to be continuous. 

Notice that the triple points A(TSG-MLC-PLC), B(OCT-TSG-PLC), and C(HCT-OCT- 
PLC) must each satisfy the conditions of two boundaries. In particular, points A and C are 
determined uniquely by dynamic properties on a homogeneous road without ramps or other 
inhomogeneities: Qf n = Q cl , Qf mp = Q c2 - Q cl , Q m = Q cv , and Q^ mp = Q out - Q cv « 

c2 Qcv 

Summarizing our results, we have presented a phase diagram of traffic states developing 
close to an on-ramp when passed by a perturbation on the freeway. Similar results are found 
for disturbances on the ramp (not displayed). Moreover, the same phase diagram is expected 
for other kinds of inhomogeneities as caused by gradients, changes in the number of lanes, 
etc. In such cases, Q Tmp corresponds to the capacity drop along the inhomogeneity. Despite 
of the complex behavior of the model, it was possible to determine the phase boundaries 
analytically and in good agreement with numerical investigations. As the analytical rela- 
tions for the phase boundaries do not contain any details of the nonlocal, gas-kinetic traffic 
model used for the simulations, they are expected to be also valid for other microscopic and 
macroscopic traffic models which have (a) the same instability diagram or, more exactly, 
critical densities p c \, p C 2, and p cv , and (b) a characteristic outflow Q ut(<5rmp, i r m P ) that is 
independent of the type of congestion at the ramp. A sensitive dependence of Q out on the 
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surrounding traffic situation could lead to different results. 

There is some evidence [§[|TB| that the homogeneous congested state can be identified 
with the observed synchronized traffic of type (i) according to the classification by Kerner 
and Rehborn M, whereas synchronized traffic of types (ii) and (iii) probably corresponds 
to the oscillating congested state. Stop-and-go waves and localized clusters have also their 
empirical counterparts f||l4| . Our findings suggest that inhomogeneities may be the generic 



reason for the formation of the different congested states |L7j . The dependence of the forming 
states on the inflows Q in and Q rmp allows to understand the observed transitions between 
localized clusters, stop-and-go waves, free and synchronized flow [j5|,|7| |i4] , |l"8l . It is also of 
importance for designing on-ramp controls. The dependence of the traffic states on the 
ramp length L rmp (because of the dependence of Q out on L rmp ) is relevant for an optimal 
dimensioning of ramps and explains why the observable traffic states are dependent on the 
observation site. 

The authors want to thank for financial support by the BMBF (research project SANDY, 
grant No. 13N7092) and by the DFG (Heisenberg scholarship He 2789/1-1). 
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FIG. 1. Phase diagram of the traffic states in the vicinity of an on-ramp as a function of 
the inflows Qi n and Q rnip on the main road and the on-ramp for fixed ramp length L = 56. The 
different states are classified at t = 135, i.e., after a sufficiently long transient period. Displayed are 
homogeneous congested traffic (O), oscillatory congested traffic (x), triggered stop-and-go traffic 
(+), moving localized clusters (□), pinned localized clusters (0), and free traffic (A). The states 
are triggered by a fully developed localized cluster travelling upstream and passing the ramp (cf. 
Fig. |3|). Solid lines indicate the theoretical phase boundaries. The dashed line represents the 
condition {Q m + (Jrmp) = Qmax that characterizes the maximum downstream flow for which a 
(possibly unstable) equilibrium solution exists. Above this line, extended congestion develops even 
without any perturbation. 
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FIG. 2. Dynamics in the flow-density space at the fixed position x = —560 upstream of the 
ramp. Dashed lines represent the flow-density diagram (i.e. the equilibrium flow-density relation), 
dotted lines the position of the critical densities p c j. (a) Shown is the transition from free traffic 
(symbols at the left side of the flow-density diagram) to homogeneous congested traffic (symbols 
at the right side). The other illustrations display (b) oscillatory congested traffic, (c) triggered 
stop-and-go traffic, and (d) a moving localized cluster. 
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FIG. 3. Spatiotemporal dynamics of typical representatives of the states depicted in Fig. |l]. 
The respective states are (a) homogeneous congested traffic (HCT), (b) oscillatory congested traffic 
(OCT), (c) triggered stop-and-go traffic (TSG), and (d) a moving localized cluster (MLC). x = 
corresponds to the middle of the ramp. 
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